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Abstract 

We introduce the notion of bilinear moment functional and study their general properties. The analogue 
of Favard's theorem for moment functionals is proven. The notion of semi-classical bilinear functionals 
is introduced as a generalization of the corresponding notion for moment functionals and motivated by 
the applications to multi-matrix random models. Integral representations of such functionals are derived 
and shown to be linearly independent. 
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1 Introduction 

The notion of moment functional is most commonly encountered as a generalization of the context of 
Orthogonal Polynomials (OP) These are generally defined as a graded polynomial orthonormal basis 
in L'^(R,d/i) where d/i is a given positive measure for which all moments 

dfj.{x)x' , (1-1) 

exist finite. The moment functional associated to such a measure is then the element C in the dual space 
of polynomials, C[x]^ defined by 

C{p{x)) / d^,p{x) , (1-2) 
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and it is uniquely characterized by its moments. The positivity of the measure impUes that we can always 
find orthogonal polynomials which are real, so that the orthogonality relation reads 



C{pra{x)pn{x)) = KSnm ■ (1-3) 

p^(x) = .t" + ©(x"-^) e R[x] , K e R^. (1-4) 
Generalizing this picture one is led to consider complex funtionals 0, i.e. whose moments are not neces- 



sarily real. The associated OPs are then defined by the same relations (1-3) where now the polynomials 
belong to the ring C[x] and ft,„ are nonzero complex numbers. 

One of the main applications of OPs is in the context of random matrices ||, ^ where they allow 
to write explicit expressions for the correlation functions of eigenvalues and of the partition function of 
these models. 

Recently ^, ^, |8| growing attention is devoted to the 2-matrix models (or the multi matrix models) 
in which the probability space is the space of couples (or n-tuples) of matrices. Also such models can 
be "solved" along lines similar to the one matrix models by finding certain bi-orthogonal polynomials 
(BOP). The probability measure is given by 



dfi{Mi,M2) = ^eT'(*^i*^^)d/ii(Mi)d//2(M2) 



(1-5) 



where Mi are N x N Hermitian matrices (usually) and the positive measures dfii are U{N) invariant. 
The relevant BOPs are then a pair of graded polynomial bases {pn{x)}, {sn(y)} "dual" to each other in 
the sense that 

dni{x)dH2iy) Pn{x)Sm{y)<3''^ = KSnm , (1-6) 

PneR[x], SneR[y] , /i„ eM^ (1-7) 



The integral in Eq. (1-6) defines a particular kind of fei- moment functional, that is an element of the dual 



to the tensor of two spaces of polynomials C[x] (8)c C[y] 

jr{p{x)\siy)) / / d^i^{x)d^l2{y)p{x)s{y)e-y , (1-8) 



provided all its bi-m,oments Hij are finite 

:=/:(x%-'') gR . (1-9) 

Generalizing this picture we now consider complex bi-momcnt functionals which are uniquely character- 
ized by their (complex) bi-moments Hij G C. 



The notion of semiclassical moment functional for a functional of the form (1-2) requires that the 
measure dfj,{x) has a density W{x) whose logaritmic derivative is a rational function of x and the support 
is a finite union of intervals. This condition can be translated into a distributional equation for the 
moment functional itself and then generalized to the complex case |l^ . 

Motivated by the applications to 2-matrix models, we are interested in the corresponding notion of 
semiclassical bi-moment functionals (which we will define properly later on) and in studying their prop- 
erties: we will produce (complex path) integral representations for them, generalizing the framework of 
^ |l3|, ^ to this situation. 
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We quickly recall that ||^, |T^, |Tl| a moment functional £ is called semi-classical if there exist two 
(minimal) fixed polynomials A{x) and B{x) with the properties that 

C {-B{x)p'{x) + A{x)p{x)) = 0, Vp(x) e C[a;] . (1-10) 

The integral representation was obtained in ||l^, |l^ |l^ : we can quickly reprove here their result (without 
details) in a different way which was not used there and which is in the line of approach of this paper. 
Consequence of the definition is that the (possibly formal) generating power series 



fc=0 



fc 



satisfies the n-th order ODE 

dz] \dz 



F{z)^0. (1-12) 



The order n is the highest of the degrees of A{x),B{x) and it is referred to -in this context- as the 
class. A distinction occurs according to the cases deg{A) < degB (Case A in Q) or deg(A) > deg(i3) 
(Case B). By looking at the recursion relation satisfied by the moments _Uk one realizes that there are 
precisely n linearly independent solutions if in Case B or n — 1 in Case A|j and hence the functionals are 



in one-to-one correspondence with the solutions of Eq. ( 1-12 ) which are analyti c at z = 0. 
It is precisely the result of jl^ that the fundamental system of solutions of Eq. ( 1-12 ) are expressible as 
Laplace integral transform of the weight density 



B{x 

(which may have also branch-points) over n distinct suitably chosen contours r^-; 

Fj{z) / dxWix)e''' . (1-14) 

In Case A one should actually reject one solution among them, i.e. the one with a singularity at the 
origin, or better consider only the linear combinations which are analytic at z = 0. 

In the present paper the bi-moment functionals we consider will rather correspond to generating 
functions in two variables satisfying an over-determined (but compatible) system of PDEs, and the 
fundamental solutions will be representable as suitably chosen double Laplace integrals. The paper is 
organized as follows: 

in Section |^ we introduce the basic objects and definitions, recalling how to explicitly construct the 
BOPs from the matrix of bimoments. We also prove that the BOPs uniquely determine the bi-moment 
functional: this is the analog in this setting of Favard's Theorem which allows to reconstruct a moment 
functional from any sequence of polynomials which satisfy a three-term recurrence relation. 
In Section ^ we introduce the definition of semiclassical functionals and then prove that (under certain 
general assumpions) they are representable as integrals of suitable 2-forms over Cartesian products of 

^In Case A and if A(x) ^ there is a linear constraint on the initial conditions for the recurrence relation, which 
decreases the dimension of solution space by one. If A{x) = then the solutions of the functional equation can be found 
easily. 
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complex paths. The starting point is the fact aheady mentioned that the generating function of bi- 
moments now depends on two variables z, w and satisfies an over-determined system of PDEs. We will 
prove the compatibility of this system (in the class of cases specified in the text) and then we will solve 
it. The solutions that we obtain (in the cases we consider) are entire functions of both variables z, w 
so that one could derive bounds on the growth of the bi-moments (the coefficients of the Taylor series 
centered ai z — — w). 

It should also be remarked that all semiclassical linear moment functionals can be recovered as a special 



case of bilinear ones (see Remark 3.1): this correspond to the fact that one-matrix models can be recovered 



from two-matrix models in which one of the measures is Gaussian. 



2 Definitions and first properties 



By bi-momcnt functional we mean a functional C on the tensor product of two copies of the space of 
polynomials 

C : C[x] (g> C[y] ^ C. (2-15) 

Although the two polynomial spaces are just copies of the same space, we use two different indeterminates 
X and y in order to distinguish them. 

Such a functional is uniquely determined by its bi-moments 



(2-16) 



It makes sense to look for bi-orthogonal polynomials. We recall their definition and some standard facts 

Definition 2.1 Two sequences of polynomials {7r„(x)},igN cind {cr„(y)}„gN of exact degree n are said to 
be biorthogonal with respect to the bi-moment functional C if 



'C(7r„|CTm) = d„ 



(2-17) 



If such two sequences exist then we denote by {pn(a;)}neN and {sn(2/)}neN the corresponding sequences 
of monic polynomials, which then satisfy 



C{pn\Sm) = hnSnm , /l„ 7^ , Vn G N. 



(2-18) 



It is an adaptation of the classical result for orthogonal polynomials to write a formula for the monic 
sequences 



Proposition 2.1 The biorthogonal polynomials exist if and only if 
A„ y^O, neN, A„ := det 



(2-19) 



Under this hypothesis the monic sequences {pn}neN and {srijneN are given by the formulas 



1 





/ Moo • 




1 


det 




■ Ml,n-1 


X 




V«,o ■ 


t^n,n—l 





(2-20) 
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Sn{y) ■■= -r— det 



Ml,0 • • • Ml,n-1 /^l,n 



V 



(2-21) 



1 



y 



y 



The proof of this simple proposition is essentially the same as for the orthogonal polynomials and it is 

left to the read er (s ee [@, |l§). 

With formula ( 2-21 ) we can also compute 

C{pn\Sm) 

The relation with the normalized polynomials is 



A„+ 



-5n 



(2-22) 



-Knix) = CnPn{x) ] cr„(y) := C„S„(2/) , (2-23) 
where the complex constants c„ and c„ are such that c„c„ 

If biorthogonal polynomials exist they in general do not satisfy a three terms recurrence relation as for 
the ordinary orthogonal polynomials: they rather satisfy recurrence relations which generally are not of 
finite bands 



-■ c„s„{y) , 

A„ + i 



XTTn{x) = 7„7r„+i(a;) + } ^ aj(n)7r„--,-(x) 

3=0 
n 



(2-24) 
(2-25) 



In the case of orthogonal polynomials the three terms recurrence relation is sufficient for reconstructing 
t he momen t functional (Favard's Theorem j^). A natural question is whether the recurrence relations 
(2-24, 2-25 ) are also sufficient for the existence of a moment bifunctional for which the two sequences are 
bi-orthogonal polynomials. Note that the specification of the numbers 7„, ai{n), i < n and 7„, Pi{n), i < n 
deter mines uniq uely the two sequences of polynomials (with the understanding that 7r_„ = = a-n) in 
Eqs. (2-24,2^2^) provided that 7„ 7^ 7^ 7„, Vn G N. The following theorem answers positively to the 
existence of the moment bifunctional 

Theorem 2.1 [Favard-like Theorem for biorthogonal polynomials] If the constants 7„,7ri do not vanish 
for all rt G N then there exists a unique moment bifunctional C for which the two sequences of polynomials 
7r„,cr„ as in Eq. (2-24, 2-25) are biorthogonal. 

Proof. As for the ordinary Favard's theorem we proceed to the construction of the bi-moments fiij = 
C{x''\y^) by induction. We introduce the associated monic polynomials by defining 



n-l 



Pn{x) 



1 .0 J. 

— 7r„(x) TT 7fc , 

Tn 

" fe=0 
_^ n — 1 

Tn J- J- 



Po{x) 

so{y) = 



(2-26) 
(2-27) 



k=0 
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The corresponding recurrence relations have the same form as in Eq. (2-24, 2-25| ) except that now the 

constants 7„ , 7„ are replaced by 1 . 

The first moment /xqo is fixed by the requirement 

1 = £(7ro|cro) = /^ooTToCTo , (2-28) 
since the polynomials tto, ctq are just nonzero constants. 

Suppose now that the moments ^ij have already been defined for i,j < N. We need then to define the 
moments ^Nj for J = 0, . . . iV — 1, and /i^jv for i = 0, . . . , — 1 and (J-nn- By imposing the orthogonality 

^ C{pn\so) ^ + ■ ■ ■ , (2-29) 

we define fJ-No, where the dots represent an expression which contains only moments already defined (i.e. 
/iio, i < N). We define by induction on j the moments ij-nj, the first having been defined above. We 
have, for j < iV — 1 

^ C{pn\sj+i) = ^iN,]+l + ■ ■ ■ , (2-30) 

where again the dots is an expression involving only previously defined moments. This defines iinj+i- 
We can repeat the arguments for the moments fJ.iN, i < N hy reversing the role of the pi's and Sj's. 
Finally the moment fi^N defined by 




^^0N \ ^ N-l 



— n ' (2-31) 



f^NO ■ ■ ■ fJ-NN 



k=0 




where the only unknown is precisely fiNN and its coefficient in the LHS does not vanish since the 
corresponding minor is just 



(2-32) 



This completes the definition of the moment bifunctional C. Q.E.D. 

We now turn our attention to some specific class of bilinear functionals C. We do not require for 
the analysis to come that the biorthogonal polynomials exist, although for applications to multimatrix 
models this is essential. In those applications the determinants A„ are proportional to the partition 
functions for the corresponding multi-matrix integrals (up to a multiplicative factor of n!) and are also 
interpretable as tau functions of KP and 2-Toda hierarchies |l8j 



3 Bilinear semiclassical functionals 

The notion of semiclassical for ordinary moment functionals and the applications to random matrices 
suggest the following 

Definition 3.1 We say that a bilinear functional C : C[x]'Sic'C[y] C is semiclassical if there exist four 
polynomials Ai{x), Bi{x) and A2{y), B2{y) of degrees ai -I- 1, 6i + 1, 02 -I- 1, 62 + 1 respectively, such that 
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the following distributional equations are fulfilled 



{D^ o Bi [x) + Ai (x)) (g) 1£ = Bi{x) (g) yC 
l®{DyoB2{y)+A2{y))C = x®B2{y)C . 



Explicitly these equations mean that, for any polynomials p{x), s{y) 

£( - Bi{x)p'{x) + Ai(a;)p(x)|s(t/)) - c(^Bi(x)p{x)\ys{y) 
C{pix)\- B2iy)s'iy) + A2{y)s{y)^ = c(^xp{x)\B2iy)siy) 



(3-33) 

(3-34) 
(3-35) 



Remark 3.1 We mentioned that any semi-classical moment functional is -in a certain sense- a special case of 
bilinear semi-classical functional. We want to clarify this relation here. Let us consider a semiclassical bifunctional 
in which A2(y) — ay and B2{y) = 1. The defining relations become 



C{~Bip' + Aip\s) = C(Bip\ys) , C{p\ - s' + ays) = C{xp\s). 



In particular for s{y) — 1 the second in Eq. (3-36) reads 



C{p\y) = -C{xp\l) 
a 



(3-36) 



(3-37) 



The claim that the reader can check directly is that the moment functional -Cr(-) := 'C(-|l) is a semiclassical 
functional in the sense explained in the introduction with A{x) — Ai{x) — ^Bi{x) and B{x) = Bi{x). It will be 
clear later on that this "reduction" corresponds to a partial integration of a Gaussian weight. 

In analogy with the orthogonal polynomials case we also define the class 

Definition 3.2 For a semi-classical bi-functional C we define its bi-class as the pair of integers 

(si, S2) = (niax(ai, 6i) + 1, max(a2, ^2) + 1) ■ (3-38) 

Note that from the definition some recurrence relations follow for the moments /iy . In order to spell 
them out we introduce the following notations for the coefficients of the polynomials Ai , Bi 



ai+l bi+1 

A,{x) = J2 «iO>' ; Bi{x) := J2 /3iO>' 
02+1 62+1 

Then the aforementioned recurrence relations are given by 

Proposition 3.1 The moments /ly of the classical bi-functional C are subject to the relations 



(3-39) 
(3-40) 



61+1 

Pl{i)l^n+J,rn+l = 


fcl + 1 

-n ^ f3l{j)fJ,n-l+j.7n 4 


ai + 1 






3=0 


62 + 1 


62+1 


02 + 1 


P2{j)thi+l,m+j = 


-m ^ (32{j)tJ.n,m-l+j 




3=0 


3=0 


3=0 



(3-41) 
(3-42) 
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Proof. 

From the definition of scmi-classicity by setting p{x) 
3^). Q.E.D. 



x" and s{y) — in the two relations (3-34 



The two recurrence relations give an overdetermined system for the moments: it is not guaranteed 
a priori that solutions exist and if they do, how many. There are now four different cases, according to 
deg(-Bi)= deg(^i); we address in the present paper the case deg(Ai) > deg(i?i), z = 1,2 (most relevant 

in the applications to random matrix models) which is the analog of Case B in and we could call 
"Case BB" . The other cases have less interesting applications in matrix models because they correspond 
to potentials (in a sense which will be clear below) which are bounded at infinity. They are certainly 
interesting from the point of view of Eqs. (3-41, 3-42| ); for example it is a simple exercise to check that if 
deg(i3i) — deg(i32) = 1 and deg(Ai) = deg(y42) = then in general no nontrivial solutions exist for Eqs 
(pft|, |3^. 

For the rest of this paper we will make the following 
Assumptions i^A) 

deg(B,) + 1 < deg(^,) , i = 1, 2. (3-43) 
Moreover in the case deg(i3i) -)- 1 = deg(Ai) and deg(i?2) + 1 = deg(yl2) we impose 



det 



ai(ai + l) /3i(5i-hl) 
/32(62 + l) a2(a2-t-l) 



^ when ai = 5i + 1, a2 = 62 + 1 



(3-44) 



Under this assumption we can prove 



Proposition 3.2 The solutions to Eqs. ( 3-41 , 3-42 ) form a vector space of dimension M := si ■ S2 = 
(ai -hi) - (02 + 1). 

Proof. The fact that the space of solutions is a vector space is obvious from the linearity of the defining 
equations. We need to prove the assertion regarding the dimension. 
We define the (possibly formal) generating function of moments 



F{z,w) 



00 ^ I. 
•s—^ z-'w 



^ jlkl 



= L 



(e^^^je^" 



(3-45) 



From the recursion relation for the moments or (equivalently) from the definition of semi-classicity, it 
follows that such function satisfies the system of PDEs 



(9^ + u;)B2(9^)- A2(9„)jF(z,u;) =0 
{d^ -f z)B^{d, )-A^{dA F{z, w) = 



(3-46) 



Conversely, any solution of this system which is analytic at z = = w provides a semi-classical bi-moment 
functional associated with the data Ai,Bi. We now count the solutions of this system. It will be clear 
later on that all the solutions are analytic at z = = w (in fact entire) so that any solution does define 
a moment functional. 

The system (|^) is a higher order overdetermined system of PDEs for the single function (or formal 
power series) F{z,w) and the compatibility is readily seen since 



(9, + w)B2{d^) ~ A2id^), (a,„ + z)Bi{d,) - Ai{d,) 



(3-47) 
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id, + w), [d^ + z)\B2{d^)Bi{d,) = (1 - l)B2{dM{dz) = 0. 



(3-48) 
(3-49) 



Now we express the system as a first order system of PDE's on tire suitable jet extension. Let us introduce 
the notation 

F^,^,{z,w):=d,^'d^''F{z,w). (3-50) 
The proof now proceeds according to the three different cases: 
Case BBl: deg(A,) > deg(B,) + 1, i = 1, 2; 

Case BB2: deg(Ai) = deg(Bi) + 1 but deg(A2) > deg(B2) + 1 (or vice-versa); 

Case BBS: deg(Ai) = deg(Bi) -f 1, deg(A2) = deg(B2) + 1. 

For convenience we set the leading coefficients of the two polynomials Ai to unity as this does not affect 
the dimension of the solution space of the system but make the formulas to come shorter to write. 
In Case BBl (a^ > &i -f 2) we can write the two first order systems for the systems 



1)2 + 1 



^ — 0, ■ ■ ■ , ai, u — . ■ ■ ■ a2 — 1 
— 0..ai — 1 



k = Q 

^2 + 1 

k = Q 



'1 + 1 



(3-51) 



J2 "2(fc)i=^„i., 



bi + l ai 

j=o _ j=0 

'i>2 + l / 

zF,_a. + 1 Y ftC^) ^Fj,k + Fj+i,k ] - Y "2(fc)-Fj,fc 

fo = V 



/i — 0, ■ ■ ■ , ai — 1, — 0, ■ ■ ■ a2 
u — 0..a2 — 1 



d,F, 



bi + l 

.1.^2 = E ^iW) 



3 = 



(3-52) 



Note that the two systems are consistent for the unknowns -F)j,,y, /i = 0, ...,ai, i> — 0, ...,02 if we have 
bi + 2 < ai, i = 1,2. 

In Case BB2 with o i = 6 1 + 1 the second system is not anymore consistent because the RHS of the third 
equation in system (3-52) contains Faj+i^o^- It must be replaced by 



^ — 0, ■ ■ ■ , ai — 1, 1/ — 0, • ■ ■ 02 
u — 0..a2 — 1 



b2 + l 



f^J,»2 + I T, l3'2{k)wFj,k - J2 "2(fc)F,,t - ^ ai(j)F,-,„2 + 
j=o L \ ;.-=o k=Q J j=0 

01-1 62 + 1 i>2 + l / "1 / 

1^ J2 &(fc)/3i(i)F, + i.fc +/3i(ai) 1^ /52(fc) ^ /3i(i)(2F,,k+F,,fc + i)-ai(j)F,„ 

i = fc = fo = \J = \ 

Finally in the Case BB3 (ai = 61 -f- 1 and 02 = ^2 + 1) we have the two systems 

9zF^,„ = F^+i,^ fj. = 0, ■ ■ ■ ,ai - 1, 1/ = 0, • ■ ■ 02 

9.F„i.„ = {0iU)(^Fj.-^ + F:.,. + l) - ai(j)F,,„) u = 0..a2 - 1 



(1 - /5i(ai)/32(Q2))a,F„i.„2 = /^iW) 



'^Fj.a^^+Y ("'ACc) - a2(fe))Fj,fc 



ai(j)F,.„, + ^ ^ /3i(i)/32(fe)a.F,,, 

j = j = fc = 



(3-53) 



(3-54) 
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and a similar system for the dw derivative. Note that in the third equation the derivatives dzFj k are 
defined by the first and second equation. 

Since now (1— /3i(ai)/32(a2)) 7^ as per the Assumption (which is (ai(ai + l)a2(a2 + 2) — /3i(ai)/32(a2)) 7^ 
if we do not assume that the polynomials Ai,A2 are monic) then the system is still well defined; on 
the other hand, if (1 — /3i(ai)/32(a2)) = then the last equation becomes a constraint^. 

It is a lengthy but straightforward check that the two systems are indeed compatible in each of the 
three cases. Since the size of the system is M = (ai + 1) • (02 + 1) = S1S2 then there are precisely M 
linearly independend solutions. Q.E.D. 



Remark 3.2 In principle we would not have to check the compatibility because we will construct later M = SiS2 
solutions to the system, which therefore will be proven to be compatible a posteriori: the point of Prop. 3.2 is 
principally that the dimension of the solution space certainly does not exceed M because that is the dimension 
of a closed system in the jet space. 



The Proposition implies that the recurrence relations ( ^-41 , 3-42| ) determine uniquely the functional C 



in terms of the moment with i = 0, . . . , ai,j = 0, . . . , 02. We need to produce AI = S1S2 linearly 
independent semiclassical functionals associated to the same data (^1,51,^2,-62) by means of integral 
representations . 

Equivalently we can produce integral representation for the M linearly independent solutions of the 
overdetermined system of PDE's ( |3-46| ) . It is precisely in this form that we will solve the problem, 
showing contextually that the generating functions are indeed entire functions of w, z. The starting point 
is to assume that such an integral representation exists: so suppose that 

F{z,w)= [ I dx^dyW{x,y)e'+y'" , (3-55) 



is a double Laplace integral representation for a solution of ( 3-46| )p[ 



Plugging such representation in the two equations in (3-46) and assuming that the contours are so chosen 
as to allow integration by parts without boundary terms, we obtain two first order equations for the 
bi- weight W{x,y) 

(Bi{x)d^ + Ai(x)+B[{x))w{x,v)^yB^{x)W{x,y) (3-56) 
(B2[y)dy + A2{y) + B'2{y))w{x,y) ^ X B2[y)W{x,y) . (3-57) 

We make the Assumption {B) that each pair [Ai^Bi) are relatively prime or at most share a factor 
{x — c) (or {y — s) ). The reason is similar to the case of ordinary semiclassical functionals. We will return 
on this genericity assumption later on. 



The two differential equations (3-56 3-57) form an overdetermined system for the biweight W{x,y) 



which is compatible and can be solved to give the only solution (up to a multiplicative nonzero constant) 

W{x,y) = Wi{x)W2{,y)e''y = exp(-yi(x) - V2{x) + xy) , (3-58) 

*We are not going to examine this case in this paper because it is more natural to study in the context of semiclassical 
functionals of type AB or AA, i.e. when deg{Ai) < deg(Bi) 

^In principle one could integrate the two-form 3/)e^^+f ^dx A Ay over any 2-cycle, but here we do not need such 

generality 
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W[{x) _ + B[{x) W!,{y) _ A^iy) + B'M 



VKi(,x) Bi(.t) ' W2{y) 

A^{x)+B[{x) 



B2{y) 



(.):=/ 
(»):=/ 



^2(?y) + i?^(y) 



(3-59) 
(3-60) 
(3-61) 



We call the two functions Vi(x), 1^2(2/) the potentials (borrowing the name from the statistical mechanic 
and random matrix context). 

A ■ -\-B' 

Note that if there are nonzero residues at the poles of g, ' then the corresponding potential have 
logarithmic singularities or poles. The general form of the biweight is 



Wi{x) := J|(x-Xj)^^ exp 



V+{x) + 



P2 



W2{y) ■■= Uiy-Y.r-e^p 



pi 

deg(Mi)<^<7i , M,{Xj)j^O 
M2{y) 



(3-62) 



fc=i 



V^+iy) 



nf=i(2/-n)'^^ 



(3-63) 



P2 



deg{M2)<J2hk , Af2(rfe)^0 



k=l 

In this formulas and in the rest of the paper Xj denote the zeroes of Bi{x), gj + 1 the corresponding 
multiplicities and —Xj are the residues at Xj of the differential dVi{x); similarly, Yh denote the zeroes of 
-62(2/), hk + l the corresponding multiplicities and —pk the residues at Yfe of the differential dV2(y)- 

The bi-class of the corresponding semiclassical bifunctional is then the total degree of the divisor of 
poles of the derivatives of the two potentials on the Riemann spheres whose affine coordinates are x and 

y 



pi 



P2 



si=di+ ^{gj + 1) , S2 = d2 + ^(/ij + 1) . 



(3-64) 



We will also use the notations Xq = 00 e P^, lo = 00 e Pj. 



3.1 The functionals 

We will define two sets of paths in the two punctured Riemann spheres P^ and Pj. We focus on the first 
sphere, the paths in the second being defined in analogous way. 

More precisely we define Si "homologically" independent paths in ¥]. \ Cx and S2 paths in P^ \ Cy where 

Cx and Cy are suitable union of cuts and points: for example the set Cx is the union of all poles and 

essential singularities of Wi {x) and cuts extending from the branchpoints to infinity. 

The reference to the homology is not in the ordinary sense: here we are considering in fact the relative 

homology of the cut-punctured sphere with prescribed sectors around the punctures. 

We first define some sectors j = 1, . . .pi, k = 0, . . . gj — 1. around the points Xj for which gj > 
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(the multiple zeroes of Bi{x)) in such a way that 

3iiVi{x)) — > +00. (3-65) 

The number of sectors for each pole is the degree of that pole in the exponential part of Wi (x) , that is 
di + I for the pole at infinity and gj for the j-th pole. Explicitly 

5, .= |x.GC, rf^ + i <arg(^)+ ^^^^ < ^^^^ I , fc = 0...di;(3-66) 

(. :e C; ' ^ + ^ < arg(. - X,) + ^^liil^liM < ^^L±i^\ , (3-67) 

A: = 0, . . . ,gj - 1, j = 1, . . . ,pi . 

These sectors are defined precisely in such a way that approaching any of the essential singularities (i.e. 
an Xj such that gj > 0) the function Wi{x) tends to zero faster than any power. 
Definition of the contours 

The definition of the contours follows directly [p^ , but we have to repeat it in both Riemann spheres. 
For the sake of completeness we recall the way they are defined. 

1. For any Xj for which there is no essential singularity (i.e. gj = 0), then we have two subcases 

(a) Corresponding to the Xj 's which are branch points or a pole {Xj E C \ N) , we take a loop 
starting at infinity in some fixed sector S^^"^ encircling the singularity and going back to infinity 
in the same sector. 

(b) For the Xj's which are regular points {Xj € N) we take a line joining Xj to infinity and 
approaching 00 in the same sector sj^^ as before. 

2. For any Xj for which there is an essential singularity (i.e. for which gj > 0) we define gj contours 
starting from Xj in the sector S^'' and returning to Xj in the next (counterclockwise) sector. 
Finally we join the singularity Xj to cxd by a path approaching 00 within the sector S^^^ chosen at 
point 1(a). 

3. For Xq := 00 we take di contours starting at Xq in tha sector S^"^ and returning at Xq in the sector 

For later convenience we also fix a sector Sl of width /3 < tt — e which contains the sector used 
above. The picture below gives an example of the typical situation, where the light grey sector represents 
Sl ■ We will make use also of the sector £ which is a sector within the dual sectoi]^ of Sl (in dark shade 
of grey in the picture): it is not difficult to realize that we can always arrange contours in such a way 
that £ is a small sector below the real positive axis (if the leading coefficient of is real and positive, 
otherwise the whole picture should be rotated appropriately). 

®Note that in our assumptions on the degrees of Ai, Bi the degrees of the essential singularity at infinity satisfy d\ > 
1 < d2 

^ We recall that for a given sector S centered around a ray a.rg{z) = qq with width A < tt, the dual sector is the 
sector centered around the ray arg(2) = tt — qq and with width tt — A. 
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We shall also require that all contours do not intersect except possibly at some Xj and that each closed 
loop should either encircle only one singularity or have one of the Xj on its support. 

The result of this procedure produces precisely si contours. By virtue of Cauchy's theorem the choice 
is largely arbitrary. 

An important feature for what follows is that when a contour Tj is closed (on the sphere IP]:), then Wi{x) 
has a singularity and/or is unbounded in the region inside Tj. We will call this property the Property 
(P)- 




Figure 1: The set of contours in the x Riemann sphere P;|,. Here we have 
three zeroes of B{x), Xi^ X2, X3, and the singularity at infinity Xq of order 
di + 1 = 5. The zero Xi has multiplicity gj + I — A and the corresponding 
essential singularity behaves like exp (x — Xi)^^, the zero X2 is a regular point 
for Wi{x), namely A2 G N and finally the zero X3 is either a branch point 
of Wi, in which case the cut extends to infinity "inside" the contour (in the 
picture), or a pole (A3 ^ N). 
We then define the fundamental functionals by 



A,(x"|2/'") 



dx A dy Wi {x)W2 (y)e^''x"y" 

J 

I = 1, . . . , si, j = 1, . . .52 ,n,TO e N. 



(3-68) 



We point out that such contours are chosen so that the corresponding functionals are defined on any 
monomials x^y'^ and such that integration by parts does not give any boundary contribution. Each such 
functional is a semi-classical functional associated to the data Ai, Bi, A2, B2 and their number is precisely 
the expected number S1S2 for the solutions of Eqs. ( |3-46 ) for the generating functions. The problem now 
is to show that they are linearly independent. 



Remark 3.3 A special care should be directed at the case di 
Indeed in this circumstance the two polynomials V-^{x) 



1, i.e. when ai 



and V2+(y) 



: 61 4- 1 and 02 = ^2 + 1. 
are just quadratic. The 
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biweight W{x,y) has then the form 



Wix, y) = exp ( - + +...)[...] . (3-69) 



The condition on the determinant (3-44) is precisely the nondegeneracy of the quadratic form — ~ ^ 



However, if \5\\cr\ < 1 then the integrals as we have defined are always divergent when two contours which stretch 
to infinity are involved. This simply means that we cannot choose the surface of integration in factorized form 

^ p{H) Y)\xt need to resort to a surface which is not factorized. 
Alternatively we can analytically continue from the region of 5, a for which the integrals are convergent. 

Some important remarks are in order. Consider the generating functions associated to these contours 

Fij{z,w):= [ dx /\dyWi{x)W2{y)e'=ye'''+y''" . (3-70) 

Jp(x)xp(!/) 

They are entire functions of z, w and hence are indeed generating functions of the bi- moment functionals 
Cij{-\-). Indeed our assumptions on the degrees guarantees that V^^ have degree at least 2, which is 
sufficient to guarantee analyticity w.r.t. z, w in the whole complex plane. 



Remark 3.4 If the index i corresponds to a bounded contour r^"^' then Fij{z,w) is a function of exponential 

' j 



type in z (similarly for w if V'f^ is bounded). 



Remark 3.5 If the index i corresponds to one of the contours F'^' defined at point 1(a) or 1(b) above, then 
Fij{z,w) is of exponential type only for z in an appropriate sector which contains the sector £ dual to the sector 
Sl. 

Before entering into the details of the proof of linear independence let us return to the Assumption (B) 
about the pairs {Ai, Bi). Suppose that -say- Ai and Bi have a common factor {x — c)^ , K > 1 and that 
they have no other common factor. That is let us suppose that 

Ai{x) = {x-cyAi{x) , Bi{x) = {x-cYBi{x) , (3-71) 
Z > < r, K -.^ mm{l, r) , 

with Ai{c) ^ Q ^ Bi{c). Then formula ( 3-59| ) would give 

T/V ^ - -^M. - (x - c)'ii + rjx - cf-^B, + {x- crB[ 
^1^"^- Wx)- {x-cYB, ' 

so that the divisor of poles od dVi{x) h as de gree less than si. Now we have two possible cases: 
(i) if ^ > r — 1 then we can recast Eq. ( |3-72| ) in the form 



Wijx) _ {x - c)'-''+iii + (r - 1)^1 + {{x - c)B,y 
Wi{x) ^ {x~c)Bi 



(3-73) 



which is equivalent to a problem in which the polynomials Ai, Bi are substituted by Ai :— (x— c)' ^^^Ai + 
{K — l)Bi and Bi := {x — c)Bi respectively, which now satisfy the assumption (F). In particular the 
definition of the contours provides the correct number of distinct contours for the new pair {Ai, Bi), that 
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is Si — r + 1 distinct contours (in the x plane). We need to recover (K — l)s2 solutions if Z > r — 1 or 
ls2 = Ks2 if Z = r — 1. 

(ii) If Z < r — 2 then we can recast Eq. ( |3-72| ) in the form 



(3-74) 



^,(^) Ai + lix - cy-^-'Bi + i^ix - cy^-'B 



Wi{x) 



{x-cy-'Bi 



now equivalent to a problem in which the polynomials Ai^Bi are substituted by Ai := Ai + K{x — 
cY^'-~^Bi and Bi :— {x — c)''^'i?i respectively, which do not have the factor {x — c) in common and 
hence satisfy the assumption (F). The definition of the contours provides the correct number of distinct 
contours for the new pair {Ai, Bi), and we need to recover Ks2 solutions. 
The next proposition shows how to recover the missing solutions. 



Proposition 3.3 If 



Ai{x) = {x - c)"^ Aiix) , Bi{x) ^ {x - c)"^ Bi{x) , K>1 



and Ai{x), Bi{x) do not vanish both at c then Eqs. (3-46) have also the solutions 



(3-75) 
(3-76) 



Proof. 

The fact that the functions ( 3-76| ) solve our system can be checked directly. 
Indeed the first eq. in (3-46) is satisfied because the differential operator reads 

(a^, + z)Bi{d,) - A,{d,) = \{d^ + z)B,{d,) ~ ii(a,)l (a, - c) 



K 



(3-77) 



and the proposed solutions are linear combination of functions of the form z'^e'^^ fr(w), r < K which are 
all in the kernel of {dz — c)^ . The second equation in (3-46) now reads 

[(9. + w)B2id^) - A2(9„)] c" / dy{y + zY c^''"+^'>W2{y) = 

= co" l^^^ dy B2{y){v + zY c^'-^+''^W2iv) + c" j ^^^^dy{^B2(y)(d, + w) - ^2(2/)^ (y + z)'e''(™+='VF2(y) = 
= e" l^^^ dy(^B2{y)(c + d,) - ^2(2/)^ (y + z)^c"<™+'^ VaCy) + c" j ^^^^dyB2(y)W2(v)(v + z)' c"" dyic""") = 
= o" l^^^ dy[B2(y)(a, + c) - ^2(2/)^ {y + zY o^^^ + ''^W2{y) + c" j ^^^^dyB2(yW2(v)(v + zY o"^ dyic""-) = 



B2{y) [d.-dy] 



(y + zY + 



{y)B2{y) - (A2{y) + B[(y))W2{y) (y + zY a 



In Case (ii) (or in Case (i) but with I — r ~ 1) these solutions are precisely the Ks2 missing solutions. 
In Case (i) with I > r only Z — 1 = i^T — 1 among the solutions (3-76) are linearly independent from those 
defined in terms of the contour integrals. To see this we write the weight 



W[{x) _ r 
Wi (x) X — c 



Ai+B[ 
B, 



(3-78) 
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Since Bi (c) 7^ then Wi (x) has a pole of order r dX x — c and can be written as 

Wi{x) = (x - cywi{x) , (3-79) 

with wi{x) analytic at a; = c and wi{c) ^ 0. The contour which comes from infinity, encircles c and goes 
back to infinity can be retracted to a circle around the pole, so that the corresponding solutions given by 
the integral representation would be 



dx Ady{x- c)-''wi{x)e'=^'+y^+'"yW2{y) = 

c| — e 

2i7r(r-l)! / dy^^-^ («;i(x)e^(^+^)) _ W2iy) 



(i) 

Such a solution is clearly an appropriate linear combination of the Ff, s j = 0,...r— 1<K— 1 with the 
nonzero coefficient wi{c) in front of pj.^ ^\ Q.E.D 



Remark 3.6 The function in Eq. (3-7(3) with j = corresponds to a moment functional £■ = Sc ^ y, where 
y is any semi-classical moment functional associated to Aijy), B2(y ) and 5c is the delta functional supported at 
a; = c on the space of polynomials C[x]. The other solutions in Eq. (3-76) with j > are also supported at c but 
are not factorized and have the form 

j 

jO^Y. ^-"^ ® ' (3-80) 

fc=0 

for suitable moment functionals 3^^. 

If there are other roots common to Ai, Bi we can repeat the procedure until we have a reduce problem 
which satisfies the Assumption (B). 

Therefore from this point on we will assume that the data {Ai, Bi, A2, B9) satisfy the Assumption 
(B). 

Theorem 3.1 The functionals Cij or -equivalently- the generating functions 

F,j{z,w):^ I dx AdyWi{x)W2{y)e'ye"'+y'^ (3-81) 

Jr<"'xr<."' 

are linearly independent 

The proof is an adaptation of with a small improvement (and a correction). We prepare a few 
lemmas. 



Lemma 3.1 [Theorem of Mergelyan (\19^, p. 367)] If E is a closed bounded set not separating the plane and if 
F{z) is continuous on E and analytic at the interior points of E, then F{z) can he uniformly approximated on E 
by polynomials. 

The next Theorem is a rephrasing of the content of ||l^ for the proof of which we refer ibidem. 

Theorem 3.2 [Miller-Shapiro Theorem] If F is a closed simple Jordan curve and F(z) is an analytic 
function (possibly with singularities and/or multivalued) in the points inside T such that the equation 

F{z)p{z)dz = (3-82) 



holds for any polynomial p(z) G (z — zo)C[z] (for some fixed zq G F), then F(z) has no singularities inside 
F and it is bounded in the interior region of and on F. 



16 



Suppose now by contradiction that there exist constants Cij not all of which zero such that 

Sl S2 

i=i j=i 

Reduction of the problem 



da; A dy T/Fi(x)W^2(2/)e='^e^^+^"' = 0. 



(3-83) 



'xr 



We claim that if Eq. (3-83) holds then we also have 



0-EE^: 



xr 



dx A Ay W^i(a;)W^2(y)e""+^™ = E E ^'A^)"^ 3^^) > 



where we have defined 



S.(z) := 



AyW2{y)& 



yw 



(3-84) 

(3-85) 
(3-86) 



(3-87) 



Indeed consider the auxiliary function of the new variable p 

Sl S2 . 

A{p-z,w):=YY^Ci3 / da;AdyW^i(x)1^2(2/)e'' 
Jr^xr'"' 

Here w play the role of parameters. This function is entire in p (because by our assumptions deg(P^j^) > 
2 and hence for all contours going to infinity the integrand goes to zero at least as exp(— |a;p — I2/P)), 
and by applying (dzdw)^ to Eq. ( ^-83 ) we have 

K 



0={d,dwr A{l;z,w)^{ — ] A{p;z,w) 



(3-88) 



p=i 



Therefore we also have A{0; z, w) = 0, Vz, w e C, which is Eq. {3-1 
This shows that proving that the functions Fij arc linearly independent is equivalent to proving that the 
two sets of functions {Si(z)}i=i.. and j{w)}j^i,,,s2 are (separately) linearly independent. 
Both the EiS and the \E'jS are now solutions of the decoupled ODEs of the same type (i.e. with linear 
coefficients) 



dz 



zBA — ]-A,{ — 



dz 



dw 



WB2 —]-A2 — 



dw 



S,(z) = 



(3-89) 
(3-90) 



Equivalently we may say that S^s and ^!jS are generating functions for the moments of semiclassical 
functionals associated to {Ai,Bi) and (^2,-62) respectively. Their linear independence was proven in 
[ p^ . Unfortunately this latter paper has a small flaw that makes one step of the proof impossible when 
deg(A,) > deg(Si) + 2 (while it is correct if deg(^,) < deg(Bi) + 2) ||2§. 

On the other side the linear independence of certain integral representation for semi-classical moment 
functionals was obtained in [ p^ ; however their definitions for the contours forces them to a procedure of 
regularization in certain cases which is elegantly bypassed by the definition of the contours in We 
prefer to fix the proof of p5|| since then we will not need any regularization. 
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3.2 Linear independence of the SjS 

In this section we prove the hnear independence of the functions . This will also prove the linear inde- 
pendence of the functions s ince they are precisely of the same form. We assume that the polynomial 



V-^{x) appearing in Eq. ( 3-62 ) has the form 

1 



This does not affect the generality of the problem inasmuch as it amounts to a rescaling of the variable x. 
To prove their linear independence we can reduce further the problem to the case where {x) = -j^x'^'^^. 
Indeed, suppose that there exist constants Aj such that 

SI n 

Wiz;vo,...,Vd):^y^Aj dxWi{x)e''^=0, (3-92) 
1=1 -^r. 



where we have emphasized the dependence on the subleading coefficients of as given in Eqs. (3-91 
3-62| ). Considering it as a function of the variables vq, ...,Vd then Eq. ( |3-92| ) implies that 



Q\a\ 



0, Va = e N"* , Vz e C. (3-93) 



Since yV(z; Sqj ^d) is clearly entire in the variables Vi, Eq. ( 3-93 ) implies that actually it does not 
depend on them. In other words if the S^s are linearly dependent with constants Ai then also the S^s 
where we "switch off" the coefhcients Vi of the potential are linearly dependent with the same constants 

A. 

Therefore it also does not affect the generality of the problem of showing linear independence to assume 
the specific form for 

We now analyze the asymptotic behavior, and we need the following definition (here given for a V^^ more 
general than the one above). 

Definition 3.3 The steepest descent contours (SDCs) for integrals of the form 

Ir{z) := / dxe-^^^^'='>+'='H{x) , (3-95) 



Jr 

with H{x) of polynomial growth at x = oo, are the contours 7^ uniquely defined, as z 00 within the 
sectors = |arg(z) e (-2(^'0)}' 

Ik := I X G C; ^iV+ix) -xz)=Q {V+{xk{z)) - zxk{z)) , 3?(V+(a;)) ^-^^ -f 00 . I , (3-96) 

where Xk{z) are the di branches of the solution to 

V^+'{x) = z , (3-97) 
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which behave as z^^i as z ^ oo in the sector, for the different determinations of the roots of z. 
Their homology class is constant as x ^ oo within the sector. 

With reference to Figure 1, the sector £ is the narrow dark-shaded dual sector of Sl (hght-shaded). 

Proposition 3.4 Let £ be the sector arg(^;) e ^— 2(^+1) 1 0^ at ^ = 00. Then the Laplace- Fourier 
transforms over the SDCs 7fc 



Fk{z) := / dxWi{x)e'=' ,k = l,...d 

J 7k 



have the following asymptotic leading behavior in the sector £ 



(3-98) 



/ \ 7-^/2^ 2A+l-d 

Fk{z) = Kd—z 2d uj 
p 

A:='^Xj , u; := e^ , 
where ii' 7^ is a constant found in the proof. 



exp 



d+1 



Z UJ 



1 + 1 - 

z 



(3-99) 
(3-100) 



Proof. 

The proof of this asymptotic is an application of the saddle point method. Writing z = \z\e^^ with the 
change x — we can rewrite the integrals 



f e-3TT-''+^+-]^(a;-X,)^^e^(^)da; = 

j=i 



\Z\<'\Z\ d 



T{x) := exp 



exp 



d+1 



Miix) 



9i 



K ^ Q ,\x\ ^ 00. 



(3-101) 



e^d^l^'^dC , (3-102) 



(3-103) 



Let us set A := and change integration variable 



s = 5(0 := 



1 



d+l 



(3-104) 



Note that the rescaling of variable leaves the contour F in the same "homology" class, so that we can 
take the contour as fixed in the ^-plane. The saddle points for this exponential are the roots of 



= S"(0 = - e*^ , 
that is the d roots of e'^. The corresponding critical values are 



(3-105) 



d+l 



wV^"^ , a;:=e2''^/^ A; = 0, . . . , d - 1. 



(3-106) 
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The map s = S{(,) is a d + 1-fold covering of the s plane by the ^-plane with square-root-type branching 
points at the sir\9). Moreover each of the d+1 sectors (around ^ = oo) for which 3?(^''+^) > is mapped 
to the single sector 

S:={sgC, -| + e<arg(s) < . (3-107) 

The inverse map £. = S,{s) is univalucd if we perform the cuts on the s plane starting at each sii'' (9) and 
going to 5i(s) = -|-oo parallel to the real axis. Such cuts are distinct for generic values of 6. We obtain 
a simply connected domain in the ,s plane (see picture). By their definiton the SDCs corresponds to 

(the two rims of) the horizontal cuts in the s-plane that go from the critical points Scr{6) to 3?(s) = -|-oo. 

The cuts are distinct if 9 (e'^^+^'''i^ ^ 9 (e'^^+'^'^^^, for j ^ k, that is away from the Stokes' lines 

at infinity 

^fe = |arg(z) = A^eizj. (3-108) 

Therefore if z approaches infinity along a ray distinct from the Stokes' lines and within the same sector 
between them, the asymptotic expansion does not change. 
Asymptotic evaluation of the steepest descent integrals 

The integrals corresponding to the steepest descent path become 



I 1^ 
\z 



e-^«e(s)Vs,kl)§ds, (3-109) 

g{s, \z\) := n (l - e^^'^'^''^^^"' , 5(^, \z\) = K 0. (3-110) 

where A := |z|~3~. The Jacobian of the change of variable has square- root types singularity at the critical 
point Sc? since the singularities (in the sense of singularity theory) of S{^) are simple and nondegenerate. 



Figure 2: The Steepest Descent contours for d = 4. The left depicts the ^-plane, the right the s-plane. 
Then the above integral becomes, upon developing the Jacobian in Puiseux series, 

\z\^ £ ""''^('' ^'^^^^'^^fs^'^ = ^^-^"^ 
I dse-'^^~^-k{s)''9{s, \z\)^==2=^={l + ...) = (3-112) 



Izl^e-'- 
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-At 



2de— "w'' 



exp 



<i+ 1 



2V7rA"' = 
.Q.E.D. 



(3-113) 
(3-114) 
(3-115) 



In particular Proposition 3.4 shows that the SDC integrals Fk are linearly independent because their 
asymptotics clearly is. 

Since the SDCs jk and the contours Tk span the same homology, we can always assume that the 5^ 
corresponding to the closed loops attached to oo are integrals over the SDC 7^ Suppose now that there 
exist constants Ai such that 

si 

Y,A,E,{z) = 0. (3-116) 
i=i 

We split the sum into two parts; the first one contains all contour integrals corresponding to the bounded 
paths, the paths joining the finite zeroes XiS to infinity, and loops attached to Xq — 00 approaching cxi 
within the sector Sl ■ We denote the subset of the corresponding indices by II ■ Now it is a simple check 
which we leave to the reader that all these integrals are of exponential type in the sector £ dual to Sj^. 
The second subset of indices Ir corresponds to the remaining contour integrals over paths whic h come 
from and return to 00 outside the sector Sl] a careful counting gives {IrI — [d/2]. The sum in ( 3-116 ) 
can be accordingly separated in 



J2 A^E,{z) = -J2 ^.2,(z) 



(3-117) 



We want to conclude that the two sides of Eq. (3-117) must vanish separately. Indeed we have remarked 
above that the LHS in (3-117) is of exponential type in the sector £. 

On the other hand we now prove that the RHS cannot be of exponential type unless each of , i ^ Ir 



vanishes. From Prop. 3.4 we deduce that among the SDC integrals there are precisely [d/2] that have a 



dominant exponential behavior of the type exp (^-j^z i aj'^j with ^{z d ui^) > in the sector £, whic 

is not of exponential type; since the SDCs can be obtained by suitable linear combinations with integer 
coefficients of the chosen contours then the [d/2\ functions S^, i £ Ir must span the same space as the 
dominant [d/2\ linearly independent SDCs in the sector £, modulo the span of S^, i E II- In formulae 



1i{Fk '■ Fk dominant in f } ~ ^{2^, Vi}modZ{S,;, i G II} ~ Z,{S,;, i e Ir/\ . 



(3-118) 



Since no nontrivial linear combination of the [d/2] dominant SDC integrals ^l-'s in £ can be of exponential 
type, the only possibility for the RHS of Eq. (3-117) to be of exponential type in the sector £ is that 



A, = 0, V^ e Ir 



Let us now focus on the terms in the LHS of Eq. (3-117). We must now prove that also Ai — 0, i £ I^. 
We can now follow pSl without hurdles. We sketch the main steps below for the sake of completeness. 



* Saying that a function is of exponential type in a given sector means that there exist constants K and C such that the 
function is bounded by Izl^e*^'^' in that sector. 
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We need to prove that 



Q{z) := ^ dxM(cc)e"" 



= ^ A^ = 0\/ieIL. 



(3-119) 



Let a be a point within the sector £ and far enough from the origin so as to leave all contours F;, i e 1^ 
to the left|^. Let us choose a contour C starting at z and going to infinity in the sector to £. Then 
we integrate Q{C)g~°''^ along C. Since e^'-^""-' Wi (x) is jointly absolutely integrable with respect to the 
arc- length on each of the Ti, i € II and C, we may interchange the order of integration to obtain 



ieiL 



Repeating this r — 1 times and then setting z = at the end, we obtain 

Y'A.f {x- ayWi{x)dx = 0, Wen. 



Let us define 



so that Eq. ( 3-121 ) is turned into 

VA,; / (x-a) ''v{x) 



v{x) := Wi{x){x-af 
Ax 



= 0, Vr e 



(3-120) 

(3-121) 
(3-122) 
(3-123) 



Let us perform the change of variable lo = ^^7^ (a homographic transformation). We denote by 7.; the 
imag es of t he contours and by f(uj) t he function v{x{lij)). 
Eq. (|3-121[ ) (or equivalently Eq. (|3-123| )) now becomes 



I dujf{uj)P{uj) = , VP e CM . 



(3-124) 



Note that in the variable uj all contours are in the finite region of the ijj-plane and the contours look like 
the ones in Figure 3 (the missing loops attached to = uj{Xo) — uj{(X)) were the contours indexed by Ij^). 




^More precisely in the half plane to the left of the perpendicular to the bi-secant of the sector £ 
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Figure 3: The contours 7^, i £ II in the uj plane. 

We denote by E the closed and bounded set in the oj plane constituted by all contours ji, i & II and 
the interiors of the closed loops. This set E satisfies the requirements of Lemma ( |3.lD . Moreover the 
contours 7^ have all the Property (p) with respect to f{uj). 

We now start proving that the AiS vanish. 
First consider a contour 7.; without interior points (i.e. those segments which join two different Xis). Let 
uj(t) be a parametric representation where t S [0, L] is the arc length parameter parameter so that uj'{t) 
is continuous and nonvanishing on [0, L]. Therefore it follows that the function 

X.M:= -'(*)' (3-125) 
is continuous on E and analytic in the interior points of E. Hence there exists a sequence of polynomials 



Pn{^) converging uniformly to Xii^) 011 ^ (by Lemma 3.1). Plugging into Eq. (3-124) and passing to 
the limit we obtain 



(3-126) 



which implies that Ai vanishes. 

Let us now consider a closed loop, say 7;. Let T(uj) be any polynomial vanishing at loq G 7/ where loq 
is the image of the (unique) zero of Bi (x) on the contour F; . Then we define 



$/(w) 



T(lu), lu £ "fi and its interior 

0, UJ E E \ {ji and its interior} 



(3-127) 



Again, 4)i{ijj) satisfies the requirement of Lemma (3.1) and hence can be approximated uniformly by a 
sequence of polynomials. Passing the limit under the integral we then obtain 



Ai I dujf{Lo)T{Lu) = 
'■yi 



VTe (cj-cjo)CM 



(3-128) 



We then use Theorem 3.2 to conclude that / should be bounded inside 7;. But this is a contradiction 
because f{uj) has the Property (p) w.r.t. 7; since v{x) — Wi{x){x — of' had the same Property w.r.t. 
the closed contour F;. This is a contradiction unless the Ai vanishes. 

Therefore we have proven that all the Ai must vanish, i.e. t he S j(z) arc linearly independent. 

Repeating for the jiw) we conclude the proof of Theorem 3.1. 



4 Conclusion 



We make a few remarks on the cases we have not considered, i.e. when deg{Ai) < deg(-Bi) for one or 
both i = 1,2. Indeed (up to some care in the definition of the contours for reasons of convergence) one 
can easily define some solutions of Eqs. (3-46) in the form of double Laplace-Fourier integrals and also 
prove their linear independence. More complicated is to produce the analog of Prop. 3.2, that is to have 
an a-priori knowledge of the dimension of the space of solutions to Eqs. (3-46): the result (which we do 
not prove here) is that there are M = S1S2 + 1 solutions. The moment recurrences (3-41 ]3^42 ) say then 
that the bifunctionals are actually M — 1 in Case AB or AI — 2 in Case AA. That is one has to give a 
criterion to select amongst the solutions to Eq. ( 3-46| ) the ones which are analytic at w = = z. We will 
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return on this point in a future publication. 

Suffices here to say that a similar problem occurs for the semi-classical moment functionals L : C[a;] C. 



As we have illustrated in the introduction the generating function satisfies Eq. 1-12, but in general not all 
solutions are analytic at z = and hence do not define any moment functional. This can be understood 
by looking at the recurrence relations satisfied by the moments: 

d k 

n^l3{j)fin+j-i = ^a(j)/x„+j , (4-129) 

where d = deg{B) > deg{A) + 1 = fc + 1. In this case the resulting d-terms recurrence relation has actually 
only d — 1 solutions because, for n = the above equation gives a constraint on the initial conditionsFj 

fe 

= ^a(j>,. (4-130) 



This should be regarded as the requirement that the solution of Eq. (1-12) be analytic at 2 = 0. 
Now, in the bilinear case we have the additional problem that the recurrence relations for the bi-moments 
are overdetermined and hence the corresponding constraint on the initial conditions must be shown to 
be compatible as well. We postpone the more detailed discussion of this problem to a future publication. 
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